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Abstract 

We review recent work on the existence and nature of cosmological singu- 
larities that can be formed during the evolution of generic as well as specific 
cosmological spacetimes in general relativity. We first discuss necessary and 
sufficient conditions for the existence of geodesically incomplete spacetimes 
based on a tensorial analysis of the geodesic equations. We then classify the 
possible singularities of isotropic globally hyperbolic universes using the Bel- 
Robinson slice energy that closely monitors the asymptotic properties of fields 
near the singularity. This classification includes all known forms of spacetime 
singularities in isotropic universes and also predicts new types. 

Key-words: Cosmological singularities, geodesic completeness, cosmologi- 
cal models 



1. Introduction 

The general issue that we address in this paper is probably the most common one 
in cosmology: 

Does the universe exist forever? 

In other words, was there a finite proper time in the past common to whole space 
before which the universe did not exist? How about in the future? Commonly, if 
yes, we say the universe has past (resp. future) singularity, otherwise the universe 
is complete. Thus completeness is the negation of singularity and vice versa. In 
order to avoid endless philosophical debates for even the simplest issues, perhaps 
the first thing one should do in a subject of this sort is to have precise mathematical 
definitions and criteria. 

The standard definition of spacetime completeness in general relativity is to say 
that spacetime is complete if it is geodesically complete that is every causal (i.e., 
timelike or null) geodesic defined in a finite subinterval can be extended to the whole 
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real line. Otherwise, spacetime is singular that is geodesically incomplete, there is 
at least one incomplete geodesic. 

The traditional criteria (cf. Hawking and Ellis (1973)) provide sufficient condi- 
tions for a spacetime to be singular, that is geodesically incomplete. These condi- 
tions basically depend on the effect of spacetime curvature on a congruence of causal 
geodesies, in particular, they predict a blow up singularity in the expansion of the 
said congruence a finite parameter value after the congruence has started to con- 
verge, that is all geodesies will meet at a point conjugate to what can be called the 
surface of last convergence. All this follows from an analysis of the equation satisfied 
by the expansion of the bunch of geodesies called the Komar- Landau- Raychaudhouri 
equation. 

A more analytic approach to the singularities of general relativity is to consider di- 
rectly the equation satisfied not by a family of geodesies as above but by each individ- 
ual geodesic, the geodesic equation. This is the approach taken in Choquet-Bruhat 
& Cotsakis (2002) wherein sufficient analytic conditions were found for geodesies to 
have infinite proper length. It is important to point out that in this approach the 
causal techniques, which are used in the traditional singularity theorems, are also 
used here, especially global hyperbolicity, however, the tensor- analytic technique is 
completely different. 

The plan of this paper is as follows. In the next Section, we review the beginnings 
of the new approach to the singularity problem and state three theorems, the first 
giving conditions for the equivalence of global hyperbolicity to slice completeness, 
another giving sufficient conditions for a spacetime to be g-complete (theorrem 3) 
and another (Theorem 2) offering a converse in the particular case of the so-called 
trivially sliced spaces. In Section 3, we state a new singularity theorem, that is 
necessary conditions for causal g-incompleteness, of particular relevance to cosmol- 
ogy which allows one to classify all possible FRW future singularities as being those 
with a non-integrable Hubble parameter. We also provide examples from the re- 
cent literature on sudden singularities and infiationary cosmology to illustrate our 
results. In Section 4, we introduce the Bel-Robinson energy as a new way towards 
a refined classification of the possible singularities in the isotropic category. This 
leads to many different types of singularities, for instance those in tahyonic cosmol- 
ogy. We provide examples which demonstrate the relevance of many of these new 
types. Lastly, we point out further consequences of our work and possible future 
extensions. 

2. Geodesies and singularities: Necessary conditions 

In this Section, we start with a very general situation. Consider a spacetime of the 
form (V, g) with V = x X, X being an interval in M (for simplicity we can take 
it to be the whole real line) and }A a smooth manifold of dimension n, in which the 
smooth, (n + 1) -dimensional, Lorentzian metric g splits as follows: 

g = -N^{ey ^gijO'e^, e^ = dt, e' = dx' + (5'dt. (i) 

Here = N(t,x'^) denotes the lapse function, j3'^{t,x^) is the shift function and the 
spatial slices Ait {= Ai x {t}) are spacelike submanifolds endowed with the time- 
dependent spatial metric gt = gijdx'^dxK We call such a spacetime a sliced space, 
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as in Cotsakis (2004). We shall assume that our sliced space {V,g) is regularly 
hyperbolic, that is the lapse function is bounded, the shift is uniformly bounded and 
the spatial metric is itself uniformly bounded below. Under these conditions we can 
prove (cf. Choquet-Bruhat k Cotsakis (2002), Cotsakis (2004)): 

Theorem 1 Let (V, g) be a regularly hyperbolic sliced space. Then the following are 
equivalent: 

1. (Alo)5'o) is a g-complete Riemannian manifold 

2. The spacetime (V, g) is globally hyperbolic 

This theorem shows that in a reasonable spacetime the basic causal notion of global 
hyperbolicity is in fact equivalent to the condition that each slice is a geodesically 
complete manifold. 

Under what conditions is global hyperbolicity equivalent to geodesic completeness 
in the original spacetime (V,5')? What is the class of sliced spaces in which such an 
equivalence holds? In a sliced space belonging to this class, in view of the results 
of the previous Theorem, geodesic completeness of the spacetime would be guessed 
simply by looking at the completeness of a slice. Let us define a trivially sliced space 
to be a spacetime with constant lapse and shift and time-independent spatial metric. 
Then we can prove the following result (Cotsakis (2004)). 

Theorem 2 Let (V, g) be a trivially sliced space. Then the following are equivalent: 

1. The spacetime {y,g) is timelike and null geodesically complete 

2. {Aio^go) is a complete Riemannian manifold 

3. The spacetime (V, g) is globally hyperbolic. 

We see that this theorem gives sufficient and necessary conditions for g-completeness 
in spacetimes of a particularly simple form, trivially sliced spaces. For more gen- 
eral situations, things can become very complicated and no general results offering 
sufficient and necessary conditions exist. 

The first result giving sufficient conditions for a spacetime to be g-complete was 
proved in Choquet-Bruhat & Cotsakis (2002). The method of proof consisted of a 
tensor analytic argument using the geodesic equations. The tangent vector m to a 
geodesic in (V, g) parametrized by arc length, or by the canonical parameter in the 
case of a null geodesic, with components dx"" / ds in the natural frame, satisfies in an 
arbitrary frame the differential equations, 

u^'Vau'^ = u'^dau'^ + ui^u'^u^ = 0. (2) 

Since = dt/ds, the zero component of the geodesic equations can be written in 
the form. 
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where we have set, 

..^f+/»< (4) 

Setting dt/ds = y and using the standard expressions for the connection coefficients, 
the geodesic equations yield 



log^= [ N-'{-doN-2diNv' + KijvV)dt, 

y{ti) Jti 



(5) 



or, using regular hyperbolicity. 



(6) 



log^ < 21ogiV-^ + f {\VN\,Nm + \K\,Nl) dt. 
Since the length (or canonical parameter extension) of the curve C is 

i It'* 

and will be infinite if ds/dt is bounded away from zero i.e., if y = dt/ds is uniformly 
bounded, we arrive at the following consequence. 

Theorem 3 // the spacetime (V, g) is globally and regularly hyperbolic and in ad- 
dition the g-norms of the space gradient of the lapse, | VA^jg, and of the extrinsic 
curvature,\K\g, are integrahle functions on the interval \t\,-\-oo), then spacetime is 
future timelike and null geodesic complete. 

Note that in this theorem, the fact that the length of the geodesies is infinite is proved 
directly and not using the equation for the expansion of the geodesic congruence as 
one does to prove the usual singularity theorems. 

3. Hubble rate and singularities 
The completeness Theorem 2 of the previous Section gives sufficient conditions for 
timelike and null geodesic completeness and therefore implies that the negations of 
each one of its conditions are necessary conditions for the existence of singularities 
(while, the singularity theorems (cf. Hawking and Ellis (1973)) provide sufficient 
conditions for this purpose). Thus this theorem is precisely what is needed for 
apphcations to the analysis of the nature of specific cosmological models known from 
other reasons to have spacetime singularities. This is the hne of thought followed in 
Cotsakis and Klaoudatou (2005). 

Since for an FRW metric, = 1 and = 0, we see that regular hyperbolicity 
is satisfied provided the scale factor a{t) is a bounded from below function of the 
proper time t on X. The condition for the integrability of the spatial gradient of 
the lapse is trivially satisfied and this is again true in more general homogeneous 
cosmologies where the lapse function N depends only on the time and not on the 
space variables. In all these cases the norm |ViV|g is zero. On the other hand, the 
condition on the norm of the extrinsic curvature is the only one which can create a 
problem. 
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For an isotropic universe \K\g' = 3(a/a)2 = SiY^, and so FRW universes in which 
the scale factor is bounded below can fail to be complete only when that norm is 
not integrable. This can happen in only one way: There is a finite time ti for which 
H fails to be integrable on the time interval [ti, oo). Since this non-integrability of 
H can be implemented in different ways, we arrive at the following result for the 
types of future singularities that can occur in isotropic universes (cf. Cotsakis and 
Klaoudatou (2005). 

Theorem 4 Necessary conditions for the existence of future singularities in globally 
hyperbolic, regularly hyperbolic FRW universes are: 

51 For each finite t, H is non-integrable on [ti,t], or 

52 H blows up in a finite time, or 

53 H is defined and integrable (that is bounded, finite) for only a finite proper time 



When does Condition 5*1 hold? It is well known that a function H{t) is integrable 
on an interval [ti,t] if H{t) is defined on [ti,t], is continuous on (ti, t) and the limits 
lim^^^+ H{t) and lim^^^- H{t) exist. Therefore there are a number of different 
ways which can lead to a singularity of the type 5*1 and such singularities are in a 
sense more subtle than the usual ones predicted by the singularity theorems. For 
instance, they may correspond to the so-called sudden, or spontaneous, singularities 
(see Barrow (2004) for this terminology) located at the right end (say ts) at which H 
is defined and finite but the left limit, lim^^^+ H{t), may fail to exist, thus making 
H non-integrable on [ti,ts], for any finite ts (which is of course arbitrary but fixed 
from the start). 

In such a model, we have a solution of the Friedmann equations for a fluid source 
with unconnected density and pressure which, in a local neighborhood of the singu- 
larity located at the time t = tg ahead, reads 



Here we take l<n<2,0<g<l,a(ts) = a<j and "^{t) is the so-called logarithmic 
psi-series which is assumed to be convergent, tending to zero as r — >■ 0. The form 
(jHl) exists as a smooth solution only on the interval (0,^^)- Also and Hg = H{ts) 
are finite but a blows up as t ^ making H continuous only on (0, ts)- In addition, 
a(0) is finite and we can extend H and define it to be finite also at 0, -f^(O) = Hq, so 
that H is defined on [0,^^]. However, since limt^o+ H{t) = ±oo, we conclude that 
this model universe implements exactly Condition 5*1 of the previous Section and 
thus H is non-integrable on [0,ts], tg arbitrary. 

This then provides an example of a singularity characterized by the fact that 
as t ts, d — > — oo while using the field equation we see that this is really a 
divergence in the pressure, p ^ oo. In particular, we cannot have in this universe a 
family of privileged observers each having an infinite proper time and finite H. A 



interval. 




(8) 
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further calculation shows that the product EapE'^'^, Eap being the Einstein tensor, is 
unbounded at ts- Hence we find that this spacetime is future geodesically incomplete. 

As another example consider a fiat FRW model, ds'^ = dt^ — a^{t)dx^, and take 
all quantities along a null geodesic with affine parameter A. Since this is conformally 
Minkowski we have d\ oc a{t)dt, or, dX = a{t)dt/a(ts), so that dX/dt = 1 for t = ts, 
where tg is a finite value of time. Borde et al (2003) then define an averaged-out 
Hubble rate in this infiating spacetime by the equation 

1 rMts) 

A[ts) - A[ti) Jx(ti) 

and so if Hav > 0, as we would expect in a truly infiationary universe, we have, 
following Cotsakis and Klaoudatou (2005), 

^ '^'^ " x{Q - x{u) A(t.) ^ ' ~ m - x{u) ^ - x{Q-x{uy 

(10) 

This shows that the affine parameter must take values only in a finite interval (other- 
wise we would get a vanishing H) which implies geodesic incompleteness. A similar 
proof is obtained for the case of a timelike geodesic. Notice that condition (fTUIl 
holds if and only if H is integrable on only a finite interval of time, thus signalling 
incompleteness according to Theorem 4. 

4. Bel-Robinson energy and the nature of singularities 
Instead of starting from the geodesic equation, we take in this Section another path, 
a tensorial in spirit approach, and consider the Bel-Robinson energy (cf. Cotsakis 
and Klaoudatou (2006)). In a sliced space this can be defined as follows. We start 
with the 2-covariant spatial electric and magnetic tensors, defined by (cf. Choquet- 
Bruhat & York (2002) 

p _ pO 

^ji = l^^'^'HikkR^^ 1 

where riijk is the volume element of the space metric g. The Bel-Robinson energy is 
then defined as the slice integral 

Bit) = \l {\E\' + \D\'+\B\' + \H\') (11) 

^ JMt 

where by |Xp = g^^ g'^^XikXji we denote the spatial norm of the 2-covariant tensor 
X. The Bel- Robinson energy is a kind of energy of the gravitational field projected 
in a sense to a slice in spacetime and as such it can be very useful in evolution 
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studies of the Einstein equations. It was in fact used in Choquet-Bruhat & Moncrief 
(2002) to prove global existence results for cosmological spacetimes. 

For an FRW universe filled with various forms of matter with metric given by 
ds^ = —dt^ + a^{t)da'^, where da^ denotes the usual time- independent metric on 
the 3-slices of constant curvature k, the norms of the magnetic parts, are 
identically zero while |£'| and \D\, the norms of the electric parts, reduce to 

\E\'^ = 3{d/af and {Dl"^ = 3 {{d/af + k/a^)\ (^2) 
Therefore the Bel-Robinson energy becomes 

B{t)^^{\E\' + \D\'), (13) 

where C is the constant volume of (or in in the case of a non-closed space) the 
3-dimensional slice at time t. 

The approach to the singularity problem introduced in Cotsakis and Klaouda- 
tou (2006) classifies the possible types of singularities that are formed in an FRW 
geometry during its cosmic evolution using the different possible combinations of 
the three main functions in the problem, namely, the scale factor a, the Hubble 
expansion rate H and the Bel Robinson energy B. If we suppose that the model has 
a finite time singularity ai, t — tg, then the possible behaviours of the functions in 
the triplet {H,a,{\E\,\D\)) in accordance with Theorem 4 are as follows: 

51 H non-integrable on [ti,t] for every t > ti 

52 H ^ 00 at ts > ti 

Ss H otherwise pathological 

Ni a^O 
A^3 a — > 00 

Bi \E\ 00, \D\ 00 
B2 \E\ < 00, \D\ 00 
Bs \E\ 00, \D\ < 00 
B4 \E\ < 00, \D\ < 00. 

The nature of a prescribed singularity is thus described completely by specifying 
the components in a triplet of the form {Si, Nj, Bi), with the indices i,j,l taking 
their respective values as above. Except some impossible cases (cf. Cotsakis and 
Klaoudatou (2006)), all other types of finite time singularities can in principle be 
formed during the evolution of FRW, matter-filled models, in general relativity or 
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other metric theories of gravity. It is interesting to note that all the standard dust 
or radiation-filled big bang singularities fall under the strongest singularity type, 
namely, the type (^i, A^i, -Bi). For example, in a flat universe flUed with dust, at 
t = we have 

a{t) oc t^^^^O, (Ni), (14) 
H oc ^ cx), (Si), (15) 
\E\'^ = 3/4H^^oo, \D\'^ ^ 3H^ ^ oo, (Bi). (16) 

Note that this scheme is organized in such a way that the character of the singulari- 
ties (i.e., the behaviour of the defining functions) becomes milder as the indices of S, 
N and B increase (this applies to both collapse and rip-type singularities). Milder 
singularities in isotropic universes are thus expected to occur as one proceeds down 
the singularity list. 

Further, we note that we can obtain a more detailed picture of the possible 
singularity formations by studying the relative asymptotic behaviours of the three 
fimctions that define the type of singularity. Using a standard notation that ex- 
presses the behaviour of two functions around the singularity at t*, we introduce the 
notion of relative strength in the singularity classification. Let /, g be two functions. 
We say that 

1. f{t) is much smaller than g{t), f{t) « g{t), if and only if 
hm,^,, f{t)/g{t) - 

2. f{t) is similar to g{t), f{t) ~ g{t), if and only if lim^^j^ fify/gif) < oo 

3. f{t) is asymptotic to g{t), f{t) g{t), if and only if lim^.^.^^ f{t)/g{t) — 1. 

Using these relations we find that for example, the standard radiation filled isotropic 
universes (with k — 0, ±1) have the asymptotic behaviours described by 

a« H « {\E\ ^ \D\), 

whereas the rest of the standard big bang singularities have 

a« H « {\E\ ~ \D\). 

This shows the precise nature of these two specific types of singularities. Many more 
examples can be found in Cotsakis and Klaoudatou (2006) (see also Cotsakis and 
Klaoudatou (2006b)). 

The Bel-Robinson energy can also be used to test g-completeness. For instance, 
consider the case of tachyonic cosmologies as in Cotsakis and Klaoudatou (2005b). 
Tachyons, phantoms, Chaplygin gases etc, represent unobserved and unknown, ten- 
sile, negative energy and/or pressure density substances, violating some or all of 
the usual energy conditions. The chief purpose of such exotic types of matter is to 
cause cosmic acceleration and drive the late phases of the evolution of the universe. 
However, they are also bound to have a nontrivial effect on the singularity structure 
of such universes. 
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A question of interest to us is under what conditions are cosmological models 
sourced by such fields g-complete? Following Cotsakis and Klaoudatou (2005b), we 
consider again the Friedman model filled with a Chaplygin gas with equation of 
state given by p = -p-"[C + - where C = A/{1 + w) - 1 and 

subject to the condition 1 + q;=1/(1 + w). The scale factor is then given by the 

form a(t) — (Cie"'^^'^ + (726*^^'^)^^^ ,r — t — to, where Ci, C2 and C3 are constants. 
Therefore we find that in the asymptotic limits r ^ and r — > 00, if tends to 
suitable constants, that is it remains finite on [to, 00) and so by our theorems the 
model is geodesically complete. We can also arrive at the same conclusion using the 
Bel- Robinson energy of this universe which in this case can be shown that it tends 
to a constant value thereby signalling completeness. 

However, this is by no means the only possible dynamical behaviour in phan- 
tom/tachyonic cosmologies. An example of g-incomplete universe is the so-called 
graduated inflationary models flrst considered in Barrow (1990). A more general 
class of models can be described by a flat space with a fluid having equation of state 
p + p — 7P^, 7 > and A < 1. This family was further analyzed in Cotsakis and 
Klaoudatou (2006) using the Bel-Robinson energy. It was found that the exact be- 
haviour of the model (which as we showed exhibits an 5*2 type singularity) described 
originally by Barrow (1990) has a more general signiflcance and the Bel-Robinson 
energy diverges making this class of models g-incomplete. 

5. Discussion 

In this paper we have reviewed the recent introduction and application of two new 
methods to tackle the singularity problem in cosmology. The first method is based 
on a direct tensorial treatment of the geodesic equations that enables one to prove 
sufficient conditions for a spacetime to be g-complcte, that is conditions under which 
geodesies have infinite length. This in conjuction with causal techniques allow one 
to show the equivalence of the deterministic concept of global hyperbolicity to the 
analytic concept of slice as well as general spacetime completeness. Also in the 
case of simple cosmological spacetimes one can prove necessary conditions for g- 
completeness, offering a converse to the completeness theorem. In this sense, a 
singular spacetime can then be tested directly as to how it becomes singular near a 
finite time singularity by checking which functions blow up and in what way at the 
finite time singularity. 

This is in sharp contrast to the traditional method of showing geodesic incom- 
pleteness in general relativity where one considers the behaviour of a geodesic con- 
gruence and indirectly shows that the geodesies in the conguence must under certain 
geometric conditions have finite length. The conditions reviewed in this paper al- 
low, in the homogeneous and isotropic complete analysis to be made for the 
finite time singularities by looking at the behaviour of the Hubble rate (extrinsic 
curvature). This leads to a classification of the relevant singularities in such models. 
In fact this approach results in an interesting demarcation and comparison of the 
singularity structures in different evolution cosmological models. 

A further advance and refinement in this approach is possible with the consider- 
ation of the Bel-Robinson energy of the spacetime in addition to the Hubble rate 
and the scale factor in such universes. In this second analytic technique the Bel- 
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Robinson energy monitors the asymptotic behaviour of the matter fields and is very 
sensitive even with shght changes in the matter content. This whole approach has 
led to a prediction of many different types of finite time singularities in cosmological 
models and their nature has become more amenable to an analytic treatment. 

There are many deep, interesting and feasible open problems that remain in 
this field and a future synthesis of causal, tensor, dynamical as well as asymptotic 
methods will prove very useful in progressing towards the goal of discovering all 
different possibilities. 
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